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ABSTRACT 

We present results of our large-area survey for z'-band dropout galaxies at z = 7 in a 1568 arcmin^ 
sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- 
and j/-band (AoS = Ifim) images with legacy data of Subaru and HST, we have identified 22 bright 
z-dropout galaxies down to y ^ 26, one of which has a spectroscopic redshift of z = 6.96 determined 
from Lya emission. The z = 7 luminosity function (LF) yields the best-fit Schechter parameters of 
0* = Omtlil X 10-3Mpc-3, M*y = -20.10 ± 0.76 mag, and a = -1.72 ± 0.65, and indicates a 
decrease from z = 6 at a > 95% confidence level. This decrease is beyond the cosmic variance in our 
two fields, which is estimated to be a factor of < 2. We have found that the cosmic star formation 
rate density drops from the peak atz = 2 — 3toz=:7 roughly by a factor of ~ 10 but not larger 
than ~ 100. A comparison with the reionization models suggests either that the Universe could not 
be totally ionized by only galaxies at z = 7, or more likely that properties of galaxies at z = 7 are 
different from those at low redshifts having, e.g., a larger escape fraction (> 0.2), a lower metallicity, 
and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, 
and the z = 6.96 galaxy with Lya emission is located at the center of an overdense region consisting 
of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized 
bubble at z = 7. 

Subject headings: galaxies: formation — galaxies: high- redshift — cosmology: observations 



1. INTRODUCTION 

Recent results from deep galaxy surveys have raised 
exciting questions about cosmic reionization and the 
early phases of galaxy formation. These surveys have 
extended the redshi ft frontie r to z ^ 6 — 10, using 
various technique s (live et al.l 120061: iStark et alJ l2007t 



lota et all 120081 : iRichard et al.l l2008t iBouwens et all 
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)09l) . These observations have reached near the epoch 
of reio nization at z > 6, which is suggested by Fan et aLl 
(|2006D who find that the Gunn-Peterson (GP) optical 
depths of SDSS QSOs significantly increase at z ~ 6. 
Although this increase can be explained by the extrap- 
olation from low z to z > 6 via adoptio n of log-nor mal 
distribution for optical depths (Becker e t al.l l2007h . it 
is also true that GP optical depths appear to increase 
significantly at around z = 6. The ionizing sources at 
z c± 6 are eflnciently ide ntified by z'-dropou t technique 
(jStiaveUi et all 12004 cf. iBunker et a"l] |2004f). However, 
the relation between reionization and ionizing sources, 
i.e. galaxies, are still not clear in our understanding 
of cosmic reionization. Since the WMAP5 polarization 
data indicate possible scenarios of an instantaneous 
reionization at z = 1 1 zfc 1 .4 and an extended reionization 
at z - 6 - 11 (jDunklev et al.ll2009f) . the ionizing photon 
production rate of galaxies at this epoch would constrain 
these models of reionization history. It is suggested 
that a substantial fraction of galaxi es have completed 
their starburst phase before z ~ 6. lEvles et all ()2007l ) 
have shown that about 40% of the IRAC detected 
z ~ 6 dropout galaxies exhibit a significant Balmer 
break, indicating that these are post-starburst galaxies 
whose major formation phase ended at z > 6. It is 



1^ Recently, ISchaerer fc de BarrosI II2009I ) argue that lEvles et al.1 
II2007I ) would overestimate the Balmer breaks by the contributions 
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important to understand when and how the progenitors 
of these post-starburst galaxies were formed. To address 
these questions, we need to study z > 7 galaxies over 
a wide range of magnitudes to derive UV luminosity 
function (LF) and density. 

Currently, only a handful of 2: ~ 7 dropout galaxy 
candidates are photometrically identified, due to obser- 
vational difficulties of identifying z ~ 7 dropout galax- 
ies wh ose detectable continuum is redshifted t o Ifim 
(e.g. iBouwens et"all l2008t lOesch et all l2009aD . Al- 
though Hubble Space Telescope (HST) images can reach 
as deep as ~ 29 AB magnitude in near-infrared (NIR) 
bands with the state-of-the-art Wide Field Camera 3 
(WFC3), the number of z ~ 7 dropout g alaxy candi- 
dates is still li mited up to ~ 10 — 20, so fa r (Oesc h et al.l 
2009b: McLure et al.ll2009bt iBunker et al.. 2009; see also 
Bouwens et al.l l2009bi r Moreover, the present studies 
cover only small areas (~ 5 arcmin^ for the HST/WFC3 
studies and ~ 100 arcmin^ even for r e cent VLT/HAWK-I 
observations; ICastellano et al.l [20091 : iHickev et al.ll2009[ l 
or small volumes b y the gravitational lensing technique 
(e.g. ~ 100 Mpc^; iStark et al.ll2007D . Ahhough grav- 
itational lensing technique can probe very faint high-z 
galaxies that can not be found in blank fields with the 
current facilities f Stark et al.' 2007: Bradlcv et al.|[200l 
[Bou wens ct al. 2009a; Zheng et al. 2009), these pencil 
beam surveys suffer from large uncertainties due to cos- 
mic variance, and miss a population of bright galaxies 
at z > 7, which plays an important role in determin- 
ing UV luminosity density for cosmic reionization and in 
understanding forming massive galaxies in their major 
star-formation phase. 

Motivated by this, we have conducted a wide-field z- 
dropout galaxy survey by performing deep z' and y-band 
imaging down to ?; = 26.0 with Subaru/Suprime-Cam 
(jMiyazaki ct al. 2002), exploiting the excellent red sensi- 
tivity of the wide-field optical camera mounted on the 
8.3m Subaru telescope. Our program realizes a very 
large volume survey for z-dropout galaxies with the rea- 
sonably deep magnitude limit, which strongly constrains 
the bright-end of UV LF at z = 7. In this paper, we 
present our bright z-dropout galaxy candidates found 
by our Subaru observations, and constrain cosmic star- 
formation history and reionization in conjunction with 
faint z-dropout galaxies identified by the recent deep 
HST/WFC3 studies. We describe our observations and 
z = 7 z-dropout galaxy samples in fj2| and f}3l respec- 
tively. We show the UV LF and UV luminosity density 
of z = 7 galaxies in 2) We discuss the cosmic star- 
formation history, cosmic reionization, and the distribu- 
tion of our z-dropout galaxies in SjSl and summarize our 
results in |JS| Throughout this paper, magnitudes are in 
the AB system. We adopt (/i, Sim, Oa, fJfc/i^, n^, erg) = 
(0.7,0.3,0.7,0.02,1.0,0.9). 

2. DATA 

2.1. Observations and Data Reduction 

We carried out deep y- and z'- band imaging sur- 
vey with Suprime-Cam in the Subaru Deep Field (SDF; 
Kashikawa erall l2004^ and GOODS-N CCiav alisco eFall 
2004al) during the last three years from 2006 to 2009. 
The y band is a custom broad-band filter centered at 

of nebular emission to their broad band photometry. 
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Fig. 1. — Response curves of y and z' bands, together with model 
spectra of a 2 ~ 7 galaxy and foreground interlopers. Black solid 
lines represent the response curves whose peaks are normalized to 
1.0, while the black dotted lines denote those showing the total sys- 
tem throughput. Dark and light gray lines are the normalized to- 
tal 3/-band response curves for the flux-weighted SDF and GOODS 
data. All of these response curves include the CCD quantum ef- 
ficiency, transmission+reflection of telescope-|-instrument optics, 
and atmospheric absorption. Red line is the model spectrum 
mimi cking the one of a typical dropout galaxy l lPapovich et al.l 
[2001D . but redshifted to z = 6.73. Green hne is a template spec- 
trum of the t ypical local elliptical (old) galaxy placed at z = 1.4 
llColeman. Wu . & Weedman 1980). Because the age of the Uni- 
verse at 2 = 1.4 is only 4 — 5 Gyr which is younger than the typical 
age of local ellipticals, the oldest galaxies at 2 = 1.4 are proba- 
bly younger than the local elliptical of this template. Thus, even 
the oldest galaxies at 2 = 1.4 should have a bluer SED than this 
elliptical template. Blue line plots a dusty-starburst spectrum at 
2 = 1.5, which is a re production of the typical dusty-starburst SED 
llCimatti et al.ir2002l) . Cyan line shows a Galactic T3 dwa rf star, 
one of the reddest Galactic dwarf stars I IKnapp et al.|[2"004l) . 

Ifim. This filter was originally made for identifying 
z ~ 6 galaxies by the two-color Lym an break selec- 
tion technique (jShimasaku et al] 120051). Although this 
filter was formerly called z^ filter in iShimasaku et aLl 
(2005), we refer to it as y filter for sake of simplicity. 
Figure [1] shows the total response of the y band fil- 
ter (central wavelength of 9860A and FWHM of 590A), 
which includes the CCD quantum efficiency, transmis- 
sion-f reflection of telescope-|-instrument optics, and at- 
mospheric absorption. Each of SDF and GOODS-N field 
is covered by one pointing of Suprimc-Cam whose FoV is 
918 arcmin^. Table [T] summarizes our observations and 
the y and z' data that we used. Before 2007, we took 
y data with MIT-Lincoln Laborator y (MIT-LL) CCDs 
installed in Suprime-Cam (Mivaza ki et al.|[2002f ). while 
2009 data were ac quired with the new fully-depleted 
Hamamatsu CCDs (jKamata et al.l [20081 ). In our anal- 
ysis, we also use archival y and z' images taken in 2001- 
2006 to make the deepest stacking of y and z' data avail- 
able in SDF and GOODS-N. The total integration time 
of our y band images are 26 and 33 hours in SDF and 
GOODS-N, respectively. The z' image of SDF is pro- 
duced via a number of variable object studies such as su- 
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pernovae (jPoznanski et al.ll2007f ). active galactic nuclei 
(Morokuma T. et al. in preparatio n), and high proper 
motion stars (|Richmond et al.ll2009[ ). 

Our data were reduced using Su prime-cam Deep 
Field REDuction p ackage fSDFRED: fYagi et alJ [200l 
lOuchi et all l2004af l. We have found that the total 
throughput in the y band is improved by nearly a fac- 
tor of 2 in 2009 with the new Hamamatsu CCDs, but 
that the shapes of y-band total response curves are al- 
most identical between the data of MIT-LL (taken before 
2007) and Hamamatsu CCDs (taken in 2009; Figure [J). 
Thus, we combine these y-band data taken with MIT-LL 
and Hamamatsu CCDs. Figure [T] displays the normalized 
total y-band response curves of the Hamamatsu CCDs 
and the flux-weighted MIT-LL-|-Hamamatsu CCDs for 
the SDF and GOODS data. The differences between all 
the response curves are negligibly small. Since we find 
a very small difference of < 0.02 mag between the re- 
sponse curves of flux-weighted SDF and GOODS-N even 
with the model spectra of extremely red objects such as 
L/T-dwarf stars and z-dropout galaxies O, we only use 
the y-band response curve of the flux-weighted SDF in 
our analysis. The 3cr sky noise of the reduced images 
are (y, z') = (26.4, 27.7) and (26.2, 26.9) magnitudes in a 
l".8-diameter circular aperture for SDF and GOODS-N, 
respectively. 

We mask contaminated areas with halos of bright stars 
and CCD blooming or low signal-to-noise (S/N) regions 
caused by dithering at around the edge of the FoV. The 
sky area after the masking is 810 and 758 arcmin^ in 
SDF and GOODS-N, respectively. The total area is 1568 
arcmin'^. The positions of y and z' are aligned based on 
hundreds of stellar objects commonly detected in both 
images. After the PSF sizes of these images are matched, 
FWHM sizes of PSFs are estimated to be ~ 0".91 and 
~ 0".87 in SDF and GOODS-N, respectively 

2.2. Matched Images 

Except for the y and z' bands, we use the legacy 
ground-based deep op tical imaging da ta for SDF 
(iKashikawa et all l200l and GOODS-N (iCaoak et all 
|200J), and the deep HST/ACS v2.0 and Spitzer/IRAC 
vO.3 available for the central :^ 160 arcmin^ field 
of GOODS-N (iGiavalisco et all [2004al : iDickinson et"an 
[200l . 15 These data are registered with stellar objects 
to match the coordinates of our Suprime-Cam y and z' 
images in SDF and GOODS-N. Table [2] summarizes all 
of imaging data used in our study. Again, the PSF sizes 
of our SDF y and z' images are matched to that of the 
pubhc SDF images with a PSF FWHM of 0".99. Note 
that the pubhc ground-based GOODS-N data include 
U , V, and i?-band images with a relatively poor seeing 
{FWHM ~ 1".5). We cannot homogenize the seeing 
sizes of GOODS-N images without a significant loss of 
S/N ratios of our y and z' images. Because our pur- 
pose of the U BVR photometry is not to measure a color 
defined by the same PSF-|-aperture but to confirm non- 
detections, we place the upper limits of detection with a 
large, 3"-diameter, aperture size for those U, V, and R 
images. 

During the observations, we took images of spectropho- 

1^ Note that HST/ACS v2.0 data are significantly deeper than 
the previously released images. 
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Fig. 2. — Colors of z' — y for various objects as a function 
of redshift. Red and blue lines represent the models of dropout 
and dusty-starburst g alaxies that are repro ductions of the typical 
SEDs shown in Papovich et al.l | |200!] ) and ICimatti et al.l (g002), 
respectivel y. Green lines indicate elliptical, Sbc, Scd, and irregular 
galaxies ([Coleman. Wu. &: W cedman 1980). Purple and cyan star 
marks are Gala ctic stars l[Gu nn fc Stryker 1983) and L1-L9/T0- 
T9 dwarf stars ( Knapp et aLfrSOOi ). Black solid line indicates our 
color selection criterion, z' — y > 1.5. The black dotted line marks 
a redshift, z = 6.5, that is roughly a lower limit of our selection. 

tomctric stand ard star of G191-B2B with y band filter 
in GOODS-N (|Okelll99nl : iBohlin et al.lll99l . The stan- 
dard star was observed 4 times under photometric con- 
dition. We calculate photometric zero-points from the 
standard star data. The photometric zero points of the 
other images, i.e. GOODS-N z', SDF y, and SDF z', are 
determined by m atching the zero po i nts with those of 
image s take n bv ICapak et al.l (I2004D . IShimasaku et al.l 
()2005f ). and IKashikawa et all ()2004l ). respectively. We 
check these photometric zero points based on colors of 
stellar objects in our field and 175 Galactic stars c alcu- 
latcd from spectra given in Gunn fc Strvkeil ([l983^. We 
find that the colors of stellar objects in o ur da ta are 
consistent with those of iGunn fc Strvke^ s (|1983f) stars 
within ~ 0.03 magnitude. 

3. CATALOGS AND SAMPLES 

3.1. Photometric Catalogs 

Source detecti on and photometr y are performed us- 
ing SExtractor (jBertin fc ArnoutsI [r996) . The y im- 
ages are chosen for our source detection. We mea- 
sure l".8-diameter aperture magnitudes, and define a 
z' — y color with these aperture magnitudes obtained 
by the dual image mode of SExtractor. We cor- 
rect the magnitudes of objects for Galactic extinction 
of E{B -V) = 0.018 (SDF) and 0.012 (GOODS-N; 
ISchlegel. Finkbeiner. fc Davis! Il998f ). A total of 63,740 
and 55,559 objects are identified down to the 4cr detec- 
tion limits in SDF (y = 26.1) and GOODS-N {y = 25.9), 
respectively. 
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3.2. Photometric Samples 

We isolate z-dropout galaxy candidates at z ^ 7 from 
foreground interlopers with our photometric data. Fig- 
ure [T] plots the model spectra of a z ~ 7 galaxy as well as 
red objects includin g low-z galaxies an d one of the red- 
dest T dwarf stars (jKnapp et al.|[200 4?) which are likely 
to be prominent interlopers in our photometric sample. 
Figure [1] shows that the spectral feature of a significant 
1216A trough for z ^ 7 galaxies is clearly different from 
the spectral shapes of foreground objects at the z'- and 
y-band wavelengths. 

In Figure [2l we present predicted z' — y colors as a 
function of redshift for these model spectra. This fig- 
ure demonstrates that no objects except z ~ 7 galaxies 
and some late-type dwarf stars have colors redder than 
z' — y > 1.5. Since the wavelength coverages of z'- and 
j/-bands are very close (Figure [J), the Lya trough of 
z ~ 7 galaxies can be more clearly distinguished from 
the red continuu m of foreground objects. Indeed, t he 
elHptical galaxy of IColeman. Wu. fc WeedmanI (|198C0 is 
bluer than z > 6.5 galaxies by A(z' — y) > 0.4 — 0.5 
even at its reddest color. A strong z' — y color cut 
can isolate z ~ 7 galaxies with a small fraction of 
foreground interlopers. On the other hand, the re- 
cent studies of z ~ 7 galaxies apply a relatively weak 
color criterion for the Lya trough, e.g. z' — Yios > 1 
or z' — Jiio > 1, with a broad band of Y/J whose 
central wavelength is redder than that of our y band 
dBouwens et a lj 120081: lOesch et a l. 2009b; Bunke r et all 
[2009i : iCastellano et al.l [20091: Iffickev et al.. ,200a) . The 



colors of z' — Fi n.-s ~ 1 and z^ — Ji i n ^ 1 fa ll below the red 
peaks of the IColeman. Wu. fc We edmari''s (19M) ellipti- 
cal galaxy (see, e.g.. Figure 2 ofOSunkcr et al. 2009), and 
these studies cannot remove such foreground red galaxies 
with a z' — Yio5 or z' — Jhq color alone. Instead, these 
studies can distinguish foreground red galaxies with their 
available de ep NIR (J, H, and /or K) photometry (e.g. 
Figure 3 of lBunker etall 120091 ). The idea of our candi- 
date selection is to discriminate foreground red galaxies 
without NIR photometry but with a color cut stronger 
than the color criterion of the other studies. In addi- 
tion to a Lya trough, a rest-frame far UV continuum be- 
low Lyman break (912A) is damped by I GM absorption 
with a large Lyman continuum opacity (jlnoue fc Iwatal 
|2008() . Because no such continuum should be identified 
at a wavelength shorter than 7000(~ 912 x [1 -I- 6.7])A 
for sources at z > 6.7, non-detection criteria should be 
given in U, B, V, and R bands whose bandpasses are 

bluer than 7000A. From the model colors and the lack of 
a far UV continuum, we define the selection criteria of 
z ~ 7 galaxies: 



y > 1.5 



U > U2a 



B > B2„ 



V > V2a 



R > 



where U2a-, B2a, V2cr, and i?2cr are the 2a limiting 
magnitudes of [/, B, V, and R images, respectively. 
The [/-band criteria is only applied to the objects in 
GOODS-N, since there are no public [/-band data in 
SDF. The 2a limiting magnitudes are {B2a , , R2a) = 
(29.4,28.6,28.8) in SDF and (C/2^, ^2^, ^2<t, -R2<t) = 
(27.3,28.0,26.7,27.0) in GOODS-N. We select z - 7 
galaxy candidates with these photometric criteria down 
to the 4ct limits of y — 26.1 and y — 25.9 in SDF and 
GOODS-N, respectively. Note that the color criterion of 
z' — y > 1.5 can be securely applied to our photometric 
catalogs down to the given y-band detection limits, be- 



cause the 2a upper limits of z' images reach 28.2 (SDF) 
and 27.4 (GOODS-N). Our candidate selection is made 
with the photometric catalogs of Suprime-Cam images 
alone both in SDF and GOODS-N. Ahhough deep HST 
images are available in the central ^ 160 arcmin'^ field of 
GOODS-N, we do not use these HST images at this stage 
to avoid making a heterogeneous sample given by the dif- 
ferent detection limits on the sky of GOODS-N. Instead, 
we take advantage of the deep HST images for confirm- 
ing non-detections of blue continuum for the candidates 
falling in the area with the HST images (see below). 

After we reject spurious sources near the spikes of 
bright sources, the residuals of sky subtraction etc. by 
visual inspection, we obtain z-dropout galaxy samples 
consisting of 15 and 7 candidates in SDF and GOODS- 
N, respectively. All of these candidates have magnitudes 
fainter than y = 25.4 and brighter than the detection 
limits in each field. Snapshots of these 22 candidates are 
displayed in Figure [3] and properties of our candidates 
are summarized in Table [3l Some of our candidates are 
detected in z'- as well as y-bands. Because the bandpass 
of z' band includes both red and blue sides of GP trough 
(1216A), the detections of faint z' counterparts are rea- 
sonable. The z'-band detections rather confirm that the 
candidates are neither spurious sources nor transients ap- 
pearing in the y images. Although we do not apply a 
criterion of non-detection in i' /I band where a UV con- 
tinuum between Lyman break (912A) and GP trough 
(1216A) falls, none of our candidates have an i'//-band 
counterpart with a flux beyond our detection limits. 

We check the spectroscopic catalo gs of SDF and 
GOOpS-N, which are o btaine d by Kashifeiwa et all 
(pOS"), IShimasaku et al.' (2003^, 'Ouchi et al." (l2004af), 
[Shimasaku et al. (2006), Kashikawa ct al. (2006, 
fv, 



'Yosh ida et al.l (1200611 live et al.l (120061). INagao et al. 
(120071 ). iHavashi et all (I2009D. a nd iLv et all i20()W loT 
SDF. and TWirth et al.l (I200I. i Reddv et all (12006). 



Barger et al.| (l2008l) . ICoh en et al.' (^2000^: ICohenI 
[199 



Steidel et all (fT996l. fT99 9. 2003), PhiUi ps et all (11997 
Lowenthal et all (jl997l ). and iDawson et all (1200 If ) 



for GOODS-N. We find that one of our z-dropout 
candidates, SDF-63544, has a spectroscopic r edshift of 
z = 6 .96, which was originally identified by live et all 
(j2006l ) in their Lya emitter (LAE) study SDF-63544 
is the first dropout galaxy at z ~ 7 with a spectro- 
scopic redshift. This object has a y-band magnitude 
of y = 25.42, which is the brightest candidate in our 
z-dropout galaxy samples. It should be noted that 
our photometric sample surely includes a real z ~ 7 
galaxy. On the other hand, we also find that none of 
ur candidates are low-z objects with a spectroscopic 
redshift. This confirms that our photometric criteria do 
not select obvious foreground objects. 

Two out of seven candidates in the GOODS-N field, 
GOODSN-152505 and GOODSN-108036, faU in the re- 
gion with the deep GOODS-N HST/ ACS and Spitzer im- 
ages. We display snapshot images in Figure [H We have 
found that neither candidate is detected in HST/ ACS 
S435, Vgoe, and 1775 bands The Spitzer images of 
GOODSN-152505 and GOODSN-108036 are confused 
by the nearby bright objects due to large PSF sizes in 



18 GOODSN-152505 is located near the edge of HST GOODS-N 
field, and not covered by tlie -B435 image. 
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Fig. 3. — Snapshots of our z-dropout candidates identified in 
SDF (left) and GOODS-N (right). Each object has images of ([/), 
B, V, R, i' , z' , and y bands with an ID number on the right. The 
size of images is 10" X 10". North is up and east is to the left. 
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Fig. 4. — Snapshots of two 2f-dropout candidates faUing in the 
area with the HST and Spitzer images in GOODS-N. From left to 
right, we display HST/ACS S435, Veoei *775i ^850, Suprime-Cam 
y, Spitzer/IRAC S.O^tm, 4.5/im, 5.8/im, S.Ofim, and Spitzer/MIPS 
24/im images. The Suprime-Cam y image is shown for comparison. 
GOODSN-152505 is not covered with the iJ435 image. The size of 
all images is 10" X 10", but that of 24^tm band is 11" X 11". North 
is up and east is to the left. 

IRAC and MIPS data. However, there are some signa- 
tures of possible counterparts for GOODSN-152505 in 
the 4.5Aini band and for GOODSN-108036 in the 3.6/im 
and 4.5yLtni bands, which would be different from the ef- 
fects of source confusion. Although photometry of these 
possible Spitzer counterparts are more or less contami- 
nated by fluxes of the nearby objects, we simply calculate 
total magnitudes from 3"-diameter aperture magnitudes 
a nd aperture correc tions for IRAC and MIPS fluxes given 
in lYan et al.l (|2005f) and the Spitzer web page respec- 
tively. Since these sources are confused by the nearby 
objects, we find that total magnitudes (ma 6,^-45) of 
GOODSN-152505 and GOODSN-108036 are fainter than 
(23.4, 24.2) and (24.8,24.4), respectively. Thus, the col- 
ors ofy — W 3 fi and y — 77x4.5 are < 1 — 2 and < 1. Because 
lEvles et all (|2007D report that their z ^ 6 galaxies with 
the IRAC detections have colors of 1 — 2 magnitudes in 
z' — 7713.6 and z' — TO4.5 (as well as J — rris.e J — 7714.5), 
the colors of our two candidates are comparable to those 
of z ~ 6 galaxies. This confirms that these two candi- 
dates are n ot extreme l y red galaxies at z ^ 2 — 3, such 
reported bv lYan et al.l ( 2004f ). whose optical to IRAC col- 



ors (z — TMa.g) exceed 3.3. Although the possible IRAC 
counterparts of our z-dropout galaxy candidates would 
indicate that there exist post-starburst galaxies even at 
z ~ 7, there remains the possibility that these two candi- 
dates are foreground interlopers. We will discuss stellar 
population of these two candidates via detailed spectral 
energy distribution model fitting after we confirm the 
redshifts of these candidates by spectroscopy. 

There are no obvious counterparts of these two candi- 
dates at longer wavelengths; IRAC 5.8/im, S.Ofim, and 
MIPS 24/im bands. The MIPS snapshot of GOODSN- 
152505 shows a source at the left side, but this MIPS 
source is a counterpart of a bright source located near 
the left corner. Because the detection limits of the three 
bands are too shallow (21 — 23 mag; Table [5]) to iden- 
tify galaxies at large distances, the non detections, again, 
confirm that these two candidates are neither extremely 
red galaxies nor AGN at low redshifts. 

4. LUMINOSITY FUNCTION 

4.1. Surface Number Densities and Detection 
Completeness 

We obtain the number counts of all y-band detected 
objects, Nsdi{m), and our z-dropout galaxy candidates, 
A'cand(777), from our photometric catalogs. We calculate 

|http://ssc. spitzer . caltech.edu /mips/apercorr / 1 
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Fig. 5. — Top panel: Detection completeness of our y-band 
images in percentage. Black and gray solid lines represent the 
completeness for l".8-diameter aperture photometry in SDF and 
GOODS-N, respectively. Bottom panel: Surface densities of our 
objects detected in our y data. Lower and upper sequences of 
points show surface densities of our 2-dropout galaxy candidates 
and all objects detected in the y images, respectively. The surface 
densities are shown with the squares (SDF) and circle (GOODS- 
N). In the upper sequence, filled black symbols are those with the 
completeness correction, while open symbols are not applied for 
the correction. In the lower sequence, filled gray symbols are sur- 
face densities that are subtracted with the numbers of contami- 
nation, and filled black symbols denote the best-estimate of our 
2-dropout galaxy surface densities with both the contamination 
and completeness corrections (see We present total magni- 

tudes from SExtractor's MAG_AUTO for the upper sequence to 
obtain the unbiased magnitude surface density relation. Although 
the photometry of 2-dropout galaxies is carried out with a 1".8- 
diameter aperture to maximize the signal-to-noise ratio of faint 
sources, the aperture corrections are applied for x-dropout galax- 
ies in the luminosity function plot of Figure [7] For the presen- 
tation purpose, we slightly shift the open and gray points along 
the abscissa. The exact magnitudes are the same as magnitudes 
of black filled points. The vertical axis on the right side indicates 
the number counts of objects, i.e. 7V/(0.5mag)/(0.2deg^), which 
approximately correspond to numbers of z-dropout galaxies iden- 
tified in each target field. Dotted line represents the best-fit power 
law to the completeness-corrected surface densities of all objects. 

the surface number densities, Sobs, by dividing Na\\{m) 
and iVcand(w) by our respective survey areas. The re- 
sults are presented in the bottom panel of FigurelH] Since 
the surface number densities of faint objects are affected 
by detection incompleteness, we estimate detection com- 
pleteness as a function of y magnitude by Monte Carlo 
simulations. We distribute 7240 artificial objects with a 
point spread function on our y-band images after adding 
photon Poisson noise, and detect them in the same man- 
ner as for the detection for our photometric catalogs with 
SExtractor. We repeat this process 20 times, and com- 
pute the ratio of recovered objects to the input objects. 
The top panel of Figure [2] shows the detection complete- 
ness of our y-band images. We find that the detection 



completeness is typically > 70% for relatively luminous 
sources with y < 25.5. The detection completeness is 
> 50% even in the faintest magnitude bins centered at 
y = 25.85 (SDF) and y = 25.65 (GOODS-N). We correct 
the surface number densities for the detection complete- 
ness, and present them in the bottom panel of Figure [5] 

4.2. Contamination 

There are four sources of contamination in our z ~ 7 
galaxy samples: 1) spurious y-band sources made of 
noise fluctuations, 2) transients, such as faint variable 
stars-|-AGN and supernovae detected in our y-band im- 
ages, 3) foreground red objects entering our samples 
due to photometric errors, 4) L/T dwarf stars satisfy- 
ing our color selection criteria. We define the numbers 
of contamination for 1), 2), 3), and 4) in our samples 
as ^cont(»7t), iVco„t(w), ^cont(w), and iV4^„t(TO),^ respec- 
tively. Below, we check the effects of contamination, and 
estimate their impacts on our z-dropout galaxy samples. 

1) Spurious sources: 

Because we push our y-band detection limits, our z ^ 7 
galaxy samples may include spurious sources made of 
noise peaks. Some of these y-band spurious sources 
could pass our selection criteria of eq. H]). To estimate 
how much spurious sources are included in our z-dropout 
galaxy samples, we carry out source detection and color 
selection same as those in ^ but with images whose 
ADU counts are multiplied by —1. We run SExtractor 
with these negative-count images, and make negative- 
y band detection catalogs. We apply the color criteria 
of eq. Ilj, and reject sources apparently made by the 
residuals of sky subtraction in the same manner as for 
the real z-dropout galaxy selection. We find and 1 
spurious z ^ 7 galaxy candidate in the negative SDF 
and GOODS-N images, respectively. We conclude that 
our samples include N^^^^{m) = (SDF) and 1.0 ± 1.0 
(GOODS-N) spurious source statistically, where the er- 
ror is estimated by Poisson statistics. 

2) Transients: 

Although our y images are the stacked data that were 
acquired in a 6-year (3-year) long period from 2003 
(2006) to 2009 for SDF (and GOODS-N; see Table[l|), we 
investigate the possibility of transients for our z-dropout 
candidates. We stack y data taken before and after 2008, 
and obtain y-band images for the two epochs. The y- 
band detection limits of the former {yi^"^^^) and the lat- 

ter (y^^'') images are (yr„r", 2/^^'') = (25.9,26.2) 
in SDF and (26.0, 25.6) in GOODS-N. We have carried 
out photometry at the positions of our z-dropout galaxy 
candidates, and found no candidates detected at the 3cr 
levels that show a significant magnitude change between 
these two epochs by > 0.5 magnitude (corresponding to 
1 — 2 sigma levels). Thus, we conclude N^^^^{m) = in 
both SDF and GOODS-N. Checking this result, we calcu- 
late an expected number of transients based on the deep 
and w i de-field transient s t udy r esults of lMorokuma et al] 
()2008f l. iMorokuma et al.l (|2008f ) present that the number 
density of transients (with the timescale greater than 
200 days) in the magnitude range of i' ~ 25.4 — 26.0 
is ~ 0.4/deg'^ Assuming the difference of magnitudes be- 
tween y and i' is negligible for transients which are mostly 
nearby objects, the expected number of transients in our 
SDF and GOODS-N samples is only - 0.1 with the de- 
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tection completeness correction. This estimate is consis- 
tent with our conclusion of no transients in our z-dropout 
galaxy samples. 

3) Foreground red objects entering our samples due to 
photometric errors: 

It is possible that some foreground objects, such as red 
galaxies at intermediate redshifts, enter our color crite- 
ria by photometric errors, although their intrinsic colors 
cannot satisfy the criterion ol z ^ y > 1.5. We make an 
input mock catalog mimicking foreground objects, and 
carry out Monte Carlo simulations with the mock catalog 
to estimate the numbers of foreground interlopers. The 
mock catalog has the same number-density distribution 
as that of all y-detected objects corrected for the detec- 
tion completeness (Filled circles at the upper sequence 
in the bottom panel of Figure [5]). In the mock catalog, 
colors of the bright (y < 24.5) objects are the same as 
those of the observed y-detected objects. Because faint 
objects have moderately large photometric errors, we do 
not use colors of objects with a magnitude fainter than 
y = 24.5. Instead, we assign color distribution of ob- 
served y-detected objects with y = 23.5 — 24.5 to the 
objects with y > 24.5 in the mock catalog. Then, we 
perform Monte-Carlo simulations in the same manner as 
i )4.11 distributing artificial objects from the mock catalog 
on our real images and detecting them with SExtractor. 
We find that expected numbers of these interlopers are 
7V3^^j(m) = 0.85 in SDF down to y = 26.1, and 0.80 in 
GOODS-N down to y = 25.9. 

4) L/T dwarf stars: 

Figure [5] indicates that none of foreground objects ex- 
cept Galactic late-type stars can satisfy our color crite- 
rion of z — y > 1.5 without photometric errors. Late-type 
stars are potentially the largest source of contamination 
in our z-dropout galaxy samples. Estimating the num- 
bers of late-type stars which contaminate our z-dropout 
samples, we carry out Monte-Carlo simulations same as 
3), but with an input mock catalog of late-type stars. 
We use the number density of L/T dwarfs as a func - 
tion of Galactic latitude presented in iRvan et al.l ()2005l ) 
who derive the number densities in 15 deep HST/ACS 
fields down to z = 26.0 at various Galactic latitude 
and fit them to the model of surface density averaged 
over Galactic longitude. The colors of L/ T dwarfs are 
calculated with the spectral templates of iKnapp et al.l 
11)04) that cover L1-L9 and T0-T9 stars. The results 
of these Monte-Carlo simulations show that the num- 
bers of late-type star contaminants are N^^^^{m) = 5.35 
down to y = 26.1 in SDF, and 2.80 down to y = 25.9 in 
GOODS-N. The uncertainties of these simulation results 
are 0.4 — 0.5 in number. 

4.3. Redshift Distribution 

We have estimated redshift distribution of our z- 
dropout galaxies, C{m,z), by Monte-Carlo simulations 
with an input mock catalog of high-z galaxies. The 
mock catalog consists of high-z galaxies whose properties 
are given with the probability distributions of i) number 
count, ii) continuum color, and iii) Lya emissivity. First, 
for the probability distribution of i), we use the surface 
number densities of our z-dropout galaxies corrected for 
contamination and completeness (black filled points in 
Figure [5]). Second, wc assume that high-z galaxies have 
the average UV continuum slope of /3 = — 2 which is 




6.5 6.7 6.9 7.1 7.3 
z 

Fig. 6. — Redshift distribution of our z-dropout galaxies in 
SDF (top) and GOODS-N (bottom). Solid lines plot the redshift 
distributions averaged over magnitudes weighted with the number 
density distributions of our z-dropout galaxies. Gray shades rep- 
resent errors of the redshift distribution estimates obtained by our 
Monte-Carlo simulations. Arrow indicates the redshift of our z- 
dropout galaxy with the spectroscopic confirmation (SDF-63544; 
Zapcc = 6.96). 

found in z 6 dropout galaxies (iStanwav et al. 120051 : 
lYan et al.|[2005HBouwens et al.ll2006ll2007ll2009d) . and 
that the probability distribution of /3 is a Gaussian func- 
tion with a standard deviation of ct/j — 0.5. We gener- 
ate spectr a with th e stellar population synthesis model 
of Bruzua l fc Chari ot (2003), and obtain galaxy spectra 
with (3 = (—3.0) — (—1.0). The model parameters are 
the same as those of the average value s estimated for 
Lyma n break galaxies (LBGs) at z = 3 (|Papovich et al.l 
|2001| ). but with a young age of 4 Myr and Calzetti dust 
extinction ranging from E{B — V) = 0.008 to 0.400, so 
that we can reproduce a set of spectra covering from the 
very blue (/3 = —3.0) to the moderately red {(i = —1.0) 
continua. We apply a GP optical depth calculated from 
iMadaul (119951 ) for the IGM absorption, and use this prob- 
ability distribution for ii). Third, for the fraction of 
Lya emitting galaxies, we assume that 30% of z = 7 
dropout galaxies have a Lya emission line with a rest- 
frame equivalent width {EWq) of > 20 A, which is the 
same as the estimates for z = 6 dropout galaxies down 
to ~ L* (|Stanway_et. al. 2p04a,b ; Vanzella^t al. 2006i; 
iDow-Hveelund et al.ll2007i : IStanwav et al.i2007i) . We add 
a Lya luminosity to these 30% of dropout galaxies with 
the EW q probability distrib ution o f LAEs at z = 5-7 de- 
rived bv lOuchi et"all (|2008f l. Since IShaplev et all (|2003l ) 
have found that about 50% of their z = 3 LBGs have 
a Lya EWq of < OA, we assume that a half of our 
z = 7 dropout galaxies have no Lya emission line. For 
the rest of 20% of high-z galaxies, we add a very weak 
EW{) = — 20 Lya line to their spectra. 

Finally, we produce an input mock catalog of high-z 
galaxies that are distributed in the radshift space ho- 
mogeneously, and apply statistical weights following the 
probability distributions of i), ii), and iii) shown above. 
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This mock catalog is used to carry out Monte-Carlo sim- 
ulations in the same manner as those in §4.21 Note that 
luminosity function around the detection limits are usu- 
ally overestimated because Poisson noise produces more 
sources with an up scattered flux than those with a down 
scattered flux at a given flux bin for objects with a typ- 
ical steep source number count (|Eddingtonl Il940f ). To 
correct this bias, we extrapolate the relation between sur- 
face number density and magnitude down to y — 27, and 
produce the mock catalog including faint sources below 
our detection limit. We obtain y-band detection catalogs 
by the Monte-Carlo simulations, and select artificial z- 
dropout galaxies with the color criteria (eq. [1]) to draw 
mock samples of z-dropout galaxies. We calculate the 
ratio of the selected objects to the input objects as a 
function of redshift, which corresponds to the redshift 
distribution of our z-dropout galaxies. Figure [5] plots 
the redshift distribution of our z-dropout galaxies, C(z), 
averaged over magnitudes with the probability distribu- 
tions of i). For both data of SDF and GOODS-N, the 
peak redshift of C{z) is z = 6.9, and 90 percent of the z- 
dropout galaxies fall in z = 6.5 — 7.1. Thus, the redshift 
window of our z-dropout selection is z = 6.9^q4. 

4.4. UV Luminosity Function 

We derive the UV luminosity function of z-dropout 
galaxies based on the numbers of our candidate galax- 
ies f §4.1|) and the contamination objects ( §4.2p . and the 
completeness as a function of redshift, i.e. redshift dis- 
tribution ( ij4.3p . We calculate the number density, n{m), 
of z-dropout galaxies at each field in a given magnitude 
bin by 

Jo -^C{m,z)dz 

where n represents the four kinds of contaminants sat- 
isfying our selection criteria in a given magnitude bin 
(n = 4; see M.2\f . and ^ is the differential cosmic vol- 
ume with an area of SDF or GOODS-N. Because we ap- 
ply a l".8-diameter aperture for z-dropout galaxy pho- 
tometry to maximize the signal-to-noise ratio, we need 
to brighten the y magnitudes with an aperture correc- 
tion to estimate total fluxes. On the other hand, our 
y-band magnitudes can be contaminated with Lya emis- 
sion lines. We should subtract the contributions of Lya 
fluxes to obtain UV-continuum magnitudes, and dim the 
y-band brightness accordingly. Moreover, we should ap- 
ply k-correction to get UV-continuum luminosities at the 
rest-frame ~ 1500A from our 1200 — 1300A luminosities 
including small portions of Lya forest below GP trough, 
which enter the bandpass of our y filter. We estimate 
a correction factor to derive total UV-continuum mag- 
nitudes at ~ 1500A from our y-band magnitudes based 
on the results of our Monte-Carlo simulations in ii4.3l 
We calculate the mean input 1500A-continuum and the 
output y-band magnitudes averaged with the statisti- 
cal weights from the probability distribution functions 
(i j4.3p . and obtain correction magnitudes of -1-0.01 and 
-f 0.05 for SDF and GOODS-N samples, respectively. We 
apply these small corrections to our UV-continuum mag- 
nitude estimates 

Figure [7] presents UV LF of z-dropout galaxies from 
our samples with the red filled squares (SDF) and cir- 
cle (GOODS) as well as our upper limits with the red 
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Fig. 7. — UV luminosity function (LF) of 2-dropout galax- 
ies, together with those at lower redshifts. Two red filled squares 
and one red filled circle present UV LF of our z-dropout galaxies 
in SDF and GOODS-N, respectively. Two red arrows with open 
circles indicate upper limits of z-dropout galaxies estimated from 
our GOODS-N data (right) and the combination of our SDF and 
GOODS-N data (left). Magenta arrow displays the lower limit 
estimated from our spectroscopically-identified z-dropout galaxy. 
Magenta inverse-triangles represent the maximal LF, i.e., the se- 
cure upper limits of our LF estimates that include no correc- 
tion for contamination. Recent measurements and upper limits 
including those from HST/WFC3 studies are also plott ed with 
r ed st ar marks (Bouwcns et al. 2008), hexagons (McLu re et al.l 

f309b') , black filled pentagons tOesch et al. 2009b), open pentagons 
Resell et al. 2009a), open diamonds ((Mannucci et al. 2007), and 
crosses [Richard et al.i j200tj ) . Although previous studies define 
Muv at the rest-frame ~ 1300 — 2000A, a k-correction between 
~ 1400A and 2000A is as small as ~ 0.07 llOesch et al.' '2009a). 
Red solid line plots our best-fit Schechter functio n of z-dropout 
galaxies at z = 7. Dotted line is UV LF at z ~ 6 IIBouwens et all 
2008). Gray shade den otes UV LF from the other various studies 
(see Ouchi et al.l l2008l and the references therein). Dashed line is 
UV LF at z ~ 3 jSteidel ct al. 1999). Ticks on the upper horizon- 
tal axis show observed broad-band magnitudes of objects at z = 7. 
At around the top of this plot, we also tick the corresponding star- 
formation rates estimated from eq. Q with no dust extinction 
correction. 

open circles. We also plot z-dropout galaxy UV L Fs de- 
rived by the deep HST NIC M0S-HWFC3 (Bouwcnsitin 
120081: lOesch et al.l |2009a,b; McLurc ct al. 2009b ) and 
the g round-based ([Richard et al. 2006 : Mannucci et al.l 
120071 ) studies, together with UV LFs at low redshifts. 
^® Because we have a spectroscopically identified galaxy 
at z = 6.96, we can also place a lower limit on the UV 
LF. We estimate a 1500A -continuum magnitude of this 
galaxy to be Myv = —21.40 ±0.31 from the y-band pho- 
tometry [y = 25.42), Lya flux (2 x 10"^'' erg s~^ cm~^; 
live et all 120061 : lOta et al.l l200l . and the mean model 
spectrum at z = 6.96 with (3 = -2.0 (S3]). We show 
this lower limit with the magenta arrow in Figure [71 

Note that our study provides, for the first time, mea- 
surements of z = 7 UV LF at the bright magnitude of 



We cannot includ e the recent results oflBunker et all (120091 ) , 
ICastellano et al.l ||2009|) . and lHickev et al.l II2009I ). because their UV 
LF measurements are not apparently presented. 
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—22 < Afuv < —21 as well as the lower limit based 
on the spectroscopically-identified z-dropout galaxy. We 
find that our UV LF at z = 7 falls significantly below 
that at z — 6. Because our z-dropout galaxy sam- 
ples are largely corrected for contamination estimated 
by the simulations ( §4.2p . these small number densi- 
ties of our z-dropout galaxies could be due to over es- 
timates of the contaminants. We derive UV LF of our 
z-dropout galaxies in the same manner as above, but 
with no contamination subtraction, and refer to these 
estimates as the maximal LF which provides conserva- 
tive upper limits. We show the maximal LF with the 
magenta inverse-triangles in Figure [T] These magenta 
inverse-triangles also fall below the z — Q UV LF mea- 
surements. Note that the z — Q UV LF measurements 
would also have the similar problems in contamination 
estimates, but our magenta inverse-triangles come below 
the gray shade area in Figure [7] that represents a vari- 
ance of z = 6 LFs derived by different studies with var- 
ious contamination estimates. Thus, we conclude with 
our LF measurements and maximal LF estimates that 
z = 7 UV LF definitely decreases from z = 6 at the 
bright end. Our conclu si on is consistent w i th the claim of 
iMannucci erall (|2007D : ICastellano et all (|2009D . More- 
over, this decreasing tendency is similar to that found at 
the faint magnitude s by the HST studi es (Mxjy > —21; 
Bouwens et"al]l2008HOesch et al.ll2009al lbl: iMcLure et alJ 



2009bl: iBunker et al ll2009ir The combination of our re- 



sults and the HST studies suggests that UV LF decreases 
at both bright and faint magnitudes from z = 6 to 7. We 
fit a Schechter function to the LF measurements from our 
and the HST studies, and obtain the best-fit parameters 
of (/)* = 0.69^0^5 X 10-3Mpc-3, M{jv = -20.10 ± 0.76 
mag, and a = — 1.72±0.65 that maximize the likelihood, 
£ = Hi p[Nohs{'mi),Ncy,p{mi\(j)* ,M^y,a)\, where p[x,^l] 
is the Gaussian distribution with a mean evaluated at 
X, and A'obs and iVoxp are, respectively, the numbers of 
galaxies within a magnitude bin of rrii from observations 
and expectations for a given set of Schechter parame- 
ters. We summarize the best-fit Schechter parameters 
in Table m and plot the best-fit function in Figure [7] with 
the red line. Although the constraints on a are weak, a 
steep s lope of a is suggestive, which is sini ilar to that at 
z < 6 (|Bouwens et al.li2008HMcLure et al.ll2009a[ ). 

Since our survey results are based on two fields, SDF 
and GOODS-N, which are well-separated on the sky, we 
can estimate the effect of the cosmic variance on the 
number densities of z-dropout galaxies. First, the num- 
bers (-|-Poisson errors) of z-dropout galaxies in samples 
of (SDF, GOODS-N) are (2 ± 1.4, 1 ± 1.0) at y < 25.5, 
(4 ± 2.0, 3 ± 1.7) at y < 25.7, and (8 ± 2.8, 7 ± 2.6) 
at y < 25.9 Since the survey areas of SDF and 
GOODS-N are comparable (~ 0.2 deg^; igU), the num- 
ber counts in these two fields are comparable within the 
errors. Second, we compare the LF measurements of 

To avoid using the dependent measurements from the same 
HST data, we exclude the results of Oesch et al. ( 2009a, b) in the 
fitting. The inclusion of Oesch et alj I12009al lh) data does not 
change our conclusions but with artificially small errors. 

Seven {= 15—8) candidates of SDF fall in a narrow magnitude 
window o{ y = 25.9 — 26.1. Because the test in i|4. 21 indicates that 
our samples include 0—1 spurious sources. Thus, the moderately 
large number of faint candidates in SDF is not due to spurious 
sources but a steep slope of number count for 2:-dropout galaxies 
at this magnitude regime. 
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Fig. 8. — Error ellipses of Schechter parameters, M* and <f>* , 
at the 1 and 2cr confidence levels. Red lines represent our results 
of 2:-dropout galaxies at z = 7. Blue, cyan, and green contours 
denote error el l ipses for galaxies at z = 4, 5, and 6 obta i ned by 
IBouwens et al.l 120081) (large contours) and IMcLure et al.l (I2009al ) 
(small contours only for z = 5 and 6). T he dotted contours ar e 
error ellipses of z = 7 galaxies estimated by IBouwens et all 120081 ). 
All of these Schechter fits are based on a ~ —1.7. Our measure- 
ments of the red contours indicate that the Schechter parameters 
of z = 7 LF differ from those of z < 6 LFs at the > 2a (i.e. > 95%) 
level. 

SDF and GOODS-N in Figure El and find that the cos- 
mic variance is smaller than a factor of ~ 2 in number 
density for the survey areas of 0.2 deg^. We estimate 
the expe cted cosmic variance w ith the analytic CDM 
model of lSheth fc TormenI ()1999[ ) from the number den- 
sity of our z-dropout galaxies down to our magnitude 
limit (0.9+o;g x 10"^ Mpc'^ at Muv < 21), assuming 
one-to-one correspondence between galaxies and dark ha- 
los. The corresponding clustering bias is 6 = 7.6 — 9.6, 
and the cosmic variance is 22-28% for an area of 0.2 deg^ 
at z = 6.5 — 7.1. We check o ur calculations with th e 
Cosmic Variance Calculator (jTrenti fc S tiaveUi 2008^). 
and find that this calculator returns a very comparable 
number of ~ 30% for the cosmic variance after subtrac- 
tion of Poisson error term. We obtain even a smaller 
cosmic variance of 9 — 12% with the CDM model of 
ISheth fc TorinCTl(fl999l) . li we do not assume the one-to- 
one correspondence, but adopt clustering bias of galaxies 
measured at a slightly lower redsh ift of z ~ 6 (6 ~ 3 — 4; 
lOuchi et all 120051 : lOverzier et al.ll2006D . Thus, the pre- 
diction of CDM model is consistent with our observa- 
tional constraints of the cosmic variance, < 2, in the two 
areas of 0.2 deg^. 

Figure El shows the error ellipses of our Schechter pa- 
rameters, M* and (jf , for a = —1.72 at the 1 and 2a con- 
fid ence levels. We also present those of LFs a t z = 6 and 
5 (Bouwens ct al. 2008; Mc Lure et ani2009a[ ) and z = 4 
fBo uwc ns et al. 2008). Note that all of these Schechter 
fits are based on a ~ —1.7 (a = [—1.66] — [—1.74] 

|http: / /solo. colorado.edu/~trenti/Cosmic Variance. html| 



10 



Ouchi et al. 



for iBouwens et al.ll2008l and a = [-1.66] - [-1.71] for 
iMcLure et al.ll2009a[ ). and that Figure [8] exclusively com- 
pares two parameters of M* and (jf . Our measurements 
(red contours in Figure[8]) indicate that the Schechter pa- 
rameters of z = 7 LF differ from those of z < 6 LFs sig- 
nificantly, and that LF decreases from z = 6 to 7 at more 
than the 2a (i.e. 95%) level. Moreover, our constraints 
oi z = 7 Schechte r parameters are consistent with those 
of IBouwens et all (|2GG8) (dotted lines in Figure [B]), but 
are stronger than those, which allow us to rule out no 
evolution at the > 95% level. Although the errors of our 
measurements are too large to distinguish between lumi- 
nosity (L*) and number (</>*) evolutions. Figure [8] imphes 
that a decrease in L* would be the dominant factor of 
the LF evolution from z = 5 — 6 to 7. 

4.5. UV Luminosity Density 

We calculate UV-luminosity densities at z = 7 from 
our UV LF with the best-fit Schechter parameters given 
in i i4.4l First, we integrate the LF down to the magni- 
tude of the faintest observed z-dropout galaxies of the 
HST/WFC3 studies, which was used for our Schechter 
parameter fitting (i.e., down to Muv = —18 or ~ O.IL*), 
and obtain the observed UV-luminosity density, Pyv- 
Because a total UV-luminosity density has to be larger 
than by the amount of the contribution from galax- 
ies fainter than the limiting magnitude, p'^^ corresponds 
to the lower limit of the UV-luminosity density. We ex- 
trapolate our LF down to i = to estimate the con- 
tribution from such faint galaxies, and obtain the UV- 
luminosity density, p^^'^'' ■ Since there exist no galaxies 
with L ~ 0, this UV-luminosity density with the LF ex- 
trapolation, Puv"''^, corresponds to the upper limit of UV 
luminosity density. We estimate p^v = 4.0^2 6 ^ 10^^ 
erg s^i Hz-iMpc-3 and p"^-^" = l.ltoi x 10^6 erg s^^ 
Hz^^Mpc^'^. Both and Puv""' are summarized in 
Table SI We discuss evolution of cosmic star-formation 
rate density f ii5.ip and reionization of the Universe fi )5.2p 
using these UV-luminosity densities. For the upper limit 
value, there are uncertainties in the faint-end slope of 
the LF, a = —1.72, which we have weakly constrained. 
If we assume the best-fit Schechter parameters but with 
a steep slope of a = —1.86, Puv''" would increase only 
by a factor of 2, which just corresponds to the Itr-upper 
error value (2.2 x lO^^erg s'^ Hz-iMpc-^) of our p'^^''\ 
On the other hand, a very steep slope of a = —1.90 could 
push up Puv''^'^ by a factor of 3, and an extreme value of 
a = -1.97 may boost ^"t'^ ^^'^ by a factor of 10. However, 
IBouwens et"all (|2008D and lMcLure et al.l (|2009aD have re- 
ported a measurements similar to ours but with a small 
uncertainty, a = -1.74 ± 0.16 and a = -1.71 ± 0.11, 
respectively, for galaxies at a slightly lower redshift of 
z = 6. Since it is unlikely that the faint-end slope evolves 
largely between z = 6 to 7, Puv"" would not be well be- 
yond its IcT-upper error value in the reasonable range of 
a. 

5. DISCUSSION 

5.1. Cosmic Star- Formation History 

We calculate cosmic star-formation rate densities 
(SFRDs) from the UV-luminosity densities, and 
Pijv"'- relation between UV lu- 

minosity and star-formation rate (SFR) given by 
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Fig. 9. — Cosmic star-formation rate density (SFRD) as a 
function of redshift. Red square and inverse-triangle represent 
the extinction corrected SFRDs integrated down to L ~ O.IL* 
(SFRD°J^,\), and to L = (SFRDcS?'"'), respectively. Open square 
and inverse-triangle are the same, but with no extinction correc- 
tion, i.e., SFRD"*"" and SFRD'^pp'^''. Magenta line with triangles 
give the allowed SFRDs at 2 = 7 which are defined by SFRD°''= 
and SFRDcorr^"^ with associated errors. We shift the magenta line 
with the triangles along the abscissa for the presentation purpose. 
Filled circles indicate total SFRDs at 2 < 6 that are compiled by 
IHopkins &: B cacom (20061). Dotted line denotes the best- fit func- 
tion t o the SFRD mea surements at 2 = — 6 with the C ole et aP 
11200 II ) formalism IHopkins fc Beacomll2003) . Open triangles show 
SFRD measurements down to 0.06L*(2 = 3) with exti nction cor- 
rection, which are recently reported by IBouwens et al.. t2009bi ). 
Because their 2 = 7 measurement is very close to our results, their 
triangle symbol at 2 = 7 is hidden behind the red filled square. 

iMadau. Pozzetti. fc DickinsonI (|1998l ): 

SFR(M0 yr-i) = Lr/y(erg s-1Hz-1)/(8 x 10^^), (3) 

where Luv is UV luminosity measured at 1500A. This re- 
lation assumes that galaxies have the Salpeter IMF with 
solar metallicity We obtain SFRD°''" = 4.95+^;^! x lO'^ 
Mq yr-i Mpc-3 and SFRD^pp"' = l-32tj;^^ x lO'^ Mq 
yr^^ Mpc^"^ from p^y and p^^^^ , respectively. 

We apply extinction correction to the SFRDs, assum- 
ing the empirical relation between the UV slope, /3, and 
extinction, ^leooi for starburst galaxies, 

^1600 = 4.43 + 1.99/? (4) 

(jMeurer. Heckman. fc Calzettil Il999( l. Following the 

observational results of z ^ 6 dr opout galaxies 

JStanwav et al. 20^ lYan et al.l 120051 : IBouwens et all 
l200a ,2007 . .2009c ). we use /3 = -2 for our z-dropout 
galaxies. We estimate the extinction-corrected SFRDs 
to be SFRD°Jj,", = 7.49t^^9^ x lO'^ M© yr"! Mpc'^ and 
SFRD^'PP^^ = 2.0012-04 X 10-2 Mq yr"! Mpc'^. Since 
SFRD"PP" is the SFRD with dust correction and extrap- 
olation of the LF down to L = 0, SFRD^p'^' is an upper 
limit of our SFRD measurements. On the other hand, 
SFRD°''*' is the SFRD estimated with neither dust ex- 
tinction correction nor LF extrapolation. Thus, SFRD"*^*^ 
is regarded as a conservative lower limit. 

Figure [9] plots the cosmic SFRDs from our measure- 
ments (squares and inverse-triangles) as well as our up- 
per and lower limits including the la errors (magenta 
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line). Figure [9] also displays the cosmic SFRDs ob- 
tained from previous studies with the assumption of 
Salpeter IMF. At z ~ — 6, w e show th e cosmic 
SFRD measurements compi led bv [ Hopkins fc Bcacom 
(jJOOe). The compilation of Hopkins & Bcacom (2006) 
covers most of SFRD measurement s made, to date, 
in va ri ous wavelength inc l uding Ha ([Glazebrook et al.l 



Il999t iTresse et aLf 120021: iHanish et a ll I2006D 

IPcrc 



infrared ( Flores et al. 1999 | : IPcrcz- Gonz alez et al.ll2005D . 



sub mm (iBarger et al. 2000 ; Hugh es et al. _ 199i 



dio ([Condon et all .2002: Ser j eant et~l |2002[). and X- 
ray ( Georgakakis et arn2003f). It also inc l udes r esults 
of iGiavalisco et al.l (l2004bD. iBunker et all (|2004D . and 
lOuchi et al.l (|2004aD for z > 4 SFRDs estimated from 
UV luminosities. Figure [9] presents the best-fit func- 
tion to the z 5, 6 SFR D measurements with th e 
ICole et al.l (|2001f ) formalism (|Hopkins fc Beacon]|[2006l ). 
We also plot SFR Ds at high redshifts re cently reported 
bv lBouwens et al.l |2008). The SFRDs of iBouwens et all 
(|2008l) are applied for extinction correction, but inte- 
grated down to only 0.2L*{z = 3). Their measure- 
ments are interpreted as lower limits of total SFRDs that 
are counterparts of our SFRD°q^j. measurement. Our 
measu rement is consistent with t hat of Bouwcn s et al.l 
In Figure m comparing the Hopk ins fc Beacoml 's 



200Q) best-fit model function (dotted line) with our con- 



straints of upper and lower limits (the magenta line with 
the error bar) , we find that the cosmic SFRD drops from 
the peak at2; = 2 — 3toz = 7 roughly by a factor of 
~ 10 (at least by a factor of > 6), but that the drop is 
not as large as a factor of ~ 100. Figure [9] indicates that 
the cosmic SFRD declines even from z = 6 to 7, which 
is originated from the decrease of UV LF from z = 6 to 
7. Note that this decline of SFRD could be weaker, if 
there exist a large populat i on of very faint galaxies, such 
suggested by Stark et all (|2007| ) for z = 9 - 10 galax- 
ies, that the present blank field surveys cannot iden- 
tify. However, the decreasing tendency of SFRD from 
z = 2 — 3 to 7 cannot be rejected, if the faint-end LF has 
a reasonable slope, a > —1.97, that changes Puv'"'^, i.e., 
SFRD^P"' by a factor of < 10 (see SH). Because this 
suppression of SFRD towards high redshifts follows the 
ev olutionary tenden c y of d ark halo mass function such 
as ISheth fc Tormenl (|1999f) . this would suggest that we 
would be witnessing the early phase of galaxy formation 
history at z = 7. 
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Fig. 10. — Emission rate of ionizing photon per comoving Mpc^ , 
on, as a function of redshift. We assume /esc — 0.2, if not oth- 
erwise specified. Square and triangle present the lower and upper 
limits of Nion at 2 ~ 7 estimated from Pu'v ^^'^ Pijy^^ > respectively. 
Three thick magenta lines with triangle/inverse-triangle represent 
the allowed Nio^ ranges at 2 = 7 for feac = 0.2, 0.05, and 1.0 
from left to right. These ranges are defined by the lower and upper 



limits of A''io 



Thin magenta lines denote the associated la errors 



corresponding to la errors of the upper and lower limit estimates. 
For the presentation purpose, we shift the red square and the ma- 
genta lines along the abscissa. The exact redshift is the same as 
the one of red triangle. Solid lines plot the model predictions of 
- ^^ion that is require d for maintaining hydrogen ionization in IGM 
lI Madau et al.l [l999[ ) with clumping factors of Chh = Ij 3, 10, 
and 30, from bottom to top. Dark gray area indicates that even 
the homogeneous Universe (Chii = 1) lacks ionizi ng photons to 
maintain hydrogen ionization of IGM in the model of lMadau et al.l 
(|T993). The photoionization rates inferred from the L ya forest are 
shown with purple circles (Bolton & Hachnelt '2007). Light gray 
shade displays constraints from Bolton & Hachnelt ( 2007) who find 
that the clumping factor is Chii ^ 3 at 2 ~ 6 (and at 2 > 6 in 
the hierarchical Universe). Diamonds and pentagons are Wion at 
2 = 4 — 6 produced by both AGN and galaxies that are calculated 
from AGN LF s (Bolton & Hachnelt 2007) and galaxy LFs given by 
lYoshida et al.l (|2006) and B ouwens et al.i (12008), respectively. Up- 
per and lower sequences of diamonds / pentagons present the upper 
and lower limits of A'ion calculated from their LFs down to the 
observation limiting luminosities and L = 0, respectively, which 
are counterparts of our and p'^^°^ measurements. Although 

our A^ion at 2 = 7 includes no AGN contribution due to no AGN 
LF measurements at 2 = 7, the contribution of AGN is probably 
smaller than < 0.08 dex at logNi^n — 50 (see text). 



5.2. Ionization Photon Budget Near the Reionization 

Epoch 

We evaluate emission rate of hydrogen ionizing pho- 
ton per comoving Mpc'^, TVion, and discuss ionizing pho- 
ton budget, i.e., whether the photon production rate of 
galaxies is larger than the recombination rate of hydro- 
gen IGM. We calculate Mon for galaxies with 



7Vio„(s-i Mpc-^) = 10 



49.7 



V1025 



V3/ VO.l 



(5) 

where is the ionizing emission density at the Ly- 
man limit in units of erg s~^Hz~^Mpc~^, as is the 
spectral index of ionizing emission, and /pan is the es- 



cape fraction of ionizing photons (iBolton fc Haehnelg 
l2007f ). We adopt = puv/6 for the spectral break 



at the Lyman limit, wh ere puy is the UV luminos- 
ity density at ~ 1500A (jMadau et al.lll999D . and as- 
su me Os = 3 that cor responds to a model spectrum 
of iLeitherer et all ()1999l ) with continuous star-formation 
history, Salpeter IMF, an d a metallicity oi Z — 0.2Zq 
(|Bolton fc Haehnelg l2007f) . We apply three /osc values; 
/e sc — 0.2 sugges t ed by the recent numerical simulations 
of lYaiimaet an (120091) ■ f., . ~ 0.05 found in LBGs at 
z ~ 3 (jShaplev et all 120061 : see also llwata et all 's l2009l 
subsequent studies), and /esc = 1-0 for the maximum 
escape fraction. Note that /esc — 0.05 is also sug- 
gested for z > 4 galaxies in the semi-analytical model 
of IWvithe et~all (|2009D . Figure [10] presents two Ni^n val- 
ues at z = 7 estimated from our p^v ^"^^ Pirv with 
/esc — 0.2. We regard these two iVion as the lower and 
upper limits (see i )4.5p . To evaluate the impacts of es- 
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cape fraction determination, we plot three magenta lines 
showing the lower and upper limits of A^ion with the three 
different escape fractions, /esc = 0.2, 0.05, and 1.0. We 
calculate A'ion at z = 4 — 6 from the best-fit Schechter 
parameters of UV LF s obta ined bv | Yo shida et al. ( 2006) , 
iBouwens et"all (|2008[ ). and lMcLure^t al.. (|2009a) . Eval- 
uating upper and lower limits of A^ion at z = 4 — 6, we 
integrate UV LFs down to L = for the upper limits 
and down to the observed magnitudes of Muv — —18 



[z = 6), Muv 



-17 (z = 5), and My y 



-16 



{z = 4) for the lower limits of 'B ouwens et al.] llOOH) and 
iMcLure et al] ()2009aD Similarly, we integrate UV LFs 
down to Muv = —20 .3 and —19.2 for the lower limits of 
lYoshida et al.l (|2006D at z = 5 and 4, respectively. These 
UV luminosity densities are used to estimate A'ion with 
eq. ([5]) and f^sc = 0.2. Th en we add ionizing photons 
from AGN given bv lBolton fc Hachnelt (2007), and plot 
them in Figure [TOl We find that iVion decreases mono- 
tonically from z = 4 to 7. 

Figure [TO] also shows Nion that is required to balance 
recombination of hy drogen IGM based on the model of 
iMadau et all (|1999f ). 



iVion(s-^ Mpc-3) = 10*^-^Chii(1 



(6) 



with clumping factors of Chii ~ 1, 3, 10, and 30. Note 
that Chii = 1 corresponds to the homogeneous Universe, 
and that the Universe at z ~ 7 should have Chii > 1 due 
to the hierarchical structure formation. If the photon 
production rates of objects fall below the model of Chii — 
1, hydrogen IGM cannot maintain the ionized state for 
any clumping factors taken at the redshift. On the other 
hand, the photoionization rates inferred from the Lya 
forest indica te that clumping factor is as small as Chii ^ 
3 at z ~ 6 ([Bolton &: HaehnelH l2007t the purple circles 
in Figure [TO]). Since Chii should monotonically decrease 
towards high redshifts in the hierarchical Universe, the 
models with Chii ^ 3 can be applied to the Universe at 
z > 6. If iVion of objects falls in or beyond the model of 
1 < Chii ^ 3 at z > 6, the ionizing photon production 
rate is high enough to maintain the ionized IGM. 

Figure [TO] presents that, in the cases of /osc — 0.2 and 
/esc = 1.0, A^ion values of z = 7 galaxies are comparable 
with those predicted by the models of 1 < Chii ^ 3. 
We do not include the AGN contribution of iVion to our 
estimate at z = 7, because no AGN UV LF data are 
available at this redshift. However, AGN contribution 
of ioni zing photon production is only logiVion — 49.3 at 
z = 6 ([Bolton fc Haehneltl [20071) . Because the comovmg 
density of luminous QSOs at z ~ 6 is 30 times smaller 
than that at z ~ 3 (jFan et al.|[200l . AGN LF prob- 
ably continues decreasing towards z ~ 7. Even if we 
assume no evolution of AGN LF from z = 6 to 7, the 
AGN contribution is negligible; only pushing 0.08 dex at 
logiVion ^ 50 in Figure [TOl 

Our results indicate that the ionizing photon budget 
just balances at z = 7 for /osc ^ 0.2. On the other 
hand, A'^ion of z = 7 galaxies with /osc — 0.05 is a fac- 
tor of three below the model of Chii = 1, which has 
the ~ 95% (2cr) confidence level. Note that /osc — 0.05 



is the measured escape fraction at z ^ 3 ([Shaplev et al.l 
l2006l ). Moreover, all of the other parameters for our A^ion 
estimates, i.e. the spectral index and break, are plausi- 
ble ones for low-z star-forming galaxies. Thus, there are 
two scenarios, a) If no properties of star-forming galax- 
ies at z = 7 are different from those at low redshifts, 
the universe could not be totally ionized by only galax- 
ies (and AGN) at z = 7 at the ~ 95% confidence level, 
b) If the properties of star-forming galaxies evolve from 
low redshifts, e.g. larger /esc (/esc ^ 0.2), lower metal- 
licity, flatter IMF, and/or less dust extinction, the Uni- 
verse at z = 7 is ionized and close to being in balance 
between the rates of ionizing photon production and re- 
combination of hydrogen IGM. If the scenario a) is true, 
the hydrogen IGM would experience a deficit of ionizing 
photon at z = 7. This implies that the Universe may 
not complete the reionization by z = 7. In this case, 
the Universe would start reioniz ation right after z = 7 
and almost complete it by z ~ 6 ([Fan et aI1l2006( ). How- 
ever, the polarization data of WMAP5 place the con- 
straints that instantaneous reionization below z = 8.2 
(= 6 .7) is rejected at the 2a (3cr) level ([Dunklev et al.l 
[2009| ). It is unlikely that reionization of the Universe 
is completed at the late epoch of z ~ 6 — 7 almost 
instantaneously. Alternatively, the scenario of a) indi- 
cates that not galaxies and AGN but some exotic pop- 
ulations or physical phenomena might play an impor- 
tant role in cosmic reionization. Although several candi- 
dates of the exotic ionizing sources are suggested, such as 
dark matte r annihilation and primord ial black holes (e.g. 
[Ricotti fc "Ostrikc r 2004i;lRicotti et al.li2008. ), no observa- 
tional evidence have been obtained. It is more likely that 
the scenario b) is correct, and that the escape fraction, 
metallicity, IMF, and dust extinction are lar ger, lower, 
flatter, and/or smaller, respectively. In fact, Inoue et al.l 
((26061) claim that /osc increases towards high redshifts at 
< z < 6 via their model calculations. Th e nurn erical 
simulations of [Razoumov fc Sommer-LarserJ ([2009( ) also 
suggest the similar evolution of /esc, which reaches 0.8 
at z = 10, and that the angular averaged escape fraction 
of /osc = 0.3 — 0.6 at z ~ 7 for galaxies similar to ours 
(SFR~ 1 — 10 Mq). Metallicity of galaxi es also impact 
on th e production of ionizing photons. [Stiavelli et al.l 
([20041 ) claim that the ionizing efficiency of a stellar pop- 
ulation increases by a factor of 3 for Salpeter IMF and 
a factor of 10 for a top-heavy IMF as th e metallicity de - 
creases from Z — Zq to Z = {) (see also ISchaerei]l2003| ) . 

Moreover, dust absorption may be important in de- 
termination o f escape fraction, as demonstrated by the 
simulations of ' Laursen et al.l ([2009 f). Because Lya emis- 
sivity would be high er at z ^ 6 than lo w-z as suggested 
by LAE studies (e.g. [Ouchi et all[2008[ ). a fiatter IMF as 
well as lower metallicity and/or less dust extinction to- 
wards high-z may be plausible. It should be noted that, 
even in the scenario b), our observational constraints are 
close to being in balance between ionizing photon pro- 
duction and recombination rates at z 7. In other 
words, we might be witnessing the final stage of reioniza- 

^'^ Since the temperature of IGM increases from 10, 000 K 
(for the solar m etallicity) to ~ 20, 000 K (for low metallicity; 



The measurements of IMcLure et al.l l|2009al ) include only 
those ai z = 6 and 5. 



lOst erbrockl 1 19891) . the low metallicity in IGM would reduce the 
recombinati on rate of IG M with the solar metallicity by a factor of 
~ 2 (iStiavelli et al.ll200l) . 
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tion with the closely balanced photon budget. It would 
provide signatu res of the n eutral fract ion evoluti o n that 
are claimed by live et alj (|2006[ ) and lOta et all (|2008D 
who find significantly less number of LAEs at 2: ~ 7 
than z ~ 6. Moreover, the scenario of b) is very consis- 
tent with the extended (z ^ 6 — 11 ) reionization picture 
suggested by iDunklev et al.l (|2009D . Note that these ar- 
guments assume that there is no emergence of a large 
population of very faint galaxies at z = 7 beyond the 
Schechter function, such claimed by IStark et al.l ()2007f ) 
for the earlie r epoc h of z = 9 — 10. On the other hand, 
ISantos et al.l ()2004l ) have found that star-formation ac- 
tivities of low-mass galaxies are suppressed at z ~ 5 via 
their Keck gravitational lensing survey. Systematic lens- 
ing surveys for z = 7 galaxies would be needed to cor- 
rectly understand the contribution from these very faint 
galaxies. Although there remain the arguments of the 
very faint galaxy population, the TVjon of z = 7 galaxies 
with /esc — 0.05 still falls below the model of Chii = 1 
with a very steep faint-end slope down to a ~ —1.90 
that would boostsp^y"''^ , i.e., the upper limit of A'ion by 
a factor of ~ 3 ( UM . 
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5.3. Distribution of Dropouts: Indication of Ionized 
Bubble? 

Figure [11] presents the sky distribution of our z- 
dropout galaxies in SDF and GOODS-N. Although the 
numbers of galaxies in each field are small, they appear 
to be clustered on the sky. The distribution in SDF 
shows possible three filamentary structures crossing at 
around the center of the field; from top to the center, 
the center to bottom right, and the center to bottom 
left. The possible filaments would extend up to ^ 60 
Mpc. Our dropout galaxies in GOODS-N are located 
mostly at the center right and top left. These sky distri- 
butions imply that our z = 7 galaxies with a bright UV 
luminosity (Muv ~ -21; SFRnodnst = 10 - 30) would 
be strongly clustered, which are similar to those at z ^ 6 
(|Ouchi et al.ll2005l : lOverzier et al.|[2006l : lLee et al.ll2009D . 
In Figure [Til we mark the positions of the UV bright- 
est galaxies with y < 25.6 among our galaxies at z = 7. 
We find that these brightest galaxies are located at the 
high density regions of z-dropout galaxies both in SDF 
and GOODS-N. In the SDF panel of Figure [IH we also 
mark the Lya emitting dropout galaxy, SDF-63544, at 
Zspec = 6.96 confirmed by spectroscopy (jive et al.ll2006l ). 
Interestingly, this Lya emitting dropout galaxy sits at 
the center of the 4 UV brightest dropout galaxies whose 
distribution exten ds by ~_30 Mpc in projection (Figure 
[TT|) . Because il)ta et al.l ()2008l ) did not confirm the other 
Lya emitter candidate found by their narrow-band sur- 
vey of z = 7 even with their deep spectroscopic data, 
this SDF-63544 may be only one with an observable Lya 
emission line at z = 7 in the SDF. It would be possible 
that an overdense region of the 4 UV brightest dropouts 
would make a well-established ionized bubble of IGM in 
the cosmic volume with a size of > 30 Mpc, and that the 
ionized bubble may allow SDF-63544 to transmit Lya to 
observers with no strong Lya damping absorption given 
by neutral hydrogen of IGM. To evaluate how much Lya 
flux is absorbed by IGM, we estimate a Lya equiva- 
lent width of the Lya emitting dropout galaxy from the 
UV continuum magnitude (—21.40 ± 0.31) and Lya flux 
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Fig. 11. — Sky distribution of our 2 = 7 dropout galaxies down to 
y = 26.1 in SDF (left) and down to y = 25.9 in GOODS-N (right). 
Red circles represent positions of all dropout galaxies, while red 
circles with a magenta circle denote bright dropouts with y < 25.6. 
Blue square indicates the position of the bright dropout with Lyo 
emission at Zspec = 6.96. Gray shades are masked areas where we 
did not use the data for our analysis. The scales on the maps are 
marked in both degrees and comoving megaparsecs in projection 
at z = 7. North is up and east is to the left in these images. 

(2 X 10^^^ erg s~^ cm~^; see ij4.4p . and obtain the rest- 
frame equivalent width of EWq = 37 ± SA. This EWq 
is about a half of the one for the case B recombination 
with no absorption of Lya (68A; Nagamine ct al. 200j; 
Ono Y. et al. 2009 submitted to MNRAS). We define 
the escape fraction of Lya emission, /^^" , by 



J CSC 



L 



Lya 
obs 



-Lya ' 



(7) 



where L^^l^ and ijjjt" are observed and intrinsic Lya 
luminosities, respectively. If we assume the case B re- 
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combination, L|fj"[ergs-i] = 1.1 x lO'*^ SFR[MQyr-i], 
and SFR from A/uv with no dust extinction correction 
(eq. ini) we estimate the Lya escape fraction of the 
Lya emitting dropout to be f^y.°' = 0.54 ± 0.12. This 
fesc" would be comparable to the one at z ~ 6 with 
no strong IGM damping wing of Lya, because the aver- 
age Lya escape fraction is /j^?" — 0. 52 — 0.54 a t z ^ 6 
estimated from the Lya opacities of iFan et aTl ()2006D , 
[Madau (1995), and Mciksin (2006|) for the case that IGM 
absorbs a blue half of symmetric Lya emission line (see 
§5.1.2 of lOuchi et aH[200l . The inferred f^/^°' would 
support the idea that the Lya emitting dropout galaxy 
sits inside a well-established ionized bubble with a neu- 
tral fraction as low as that at z ~ 6. 

If this naive physical picture is correct, the size of ion- 
ized bubble would be larger than 30 Mpc at z = 7. 
Assuming that this is a lower limit of the character- 
istic bubble size , we find that the analytic models of 
iFurlanetto et al.l (j2006) would suggest an upper limit of 
neutral fraction of xm ^ 20% at z = 7. This small up- 
per limit of xm may indicate that the Universe is not 
fully neutral at z = 7. The combination of reionization 
models and LAE LF (-|-clustering) gives constra ints on 
neutral fraction of x m 5, 50% a t z = 6 .5 (.Dii kstra et al.l 



I2007t iMcQuinn et"a l. 2007; Ilic v et all 
10 — 60% at z = 7 (.Kobavashi et al 
120081 ). Our possible upper limit of xhi ^ 20% at z 



20081) and xm < 
1mm lota et all 
7 



is consistent with those of previous results from the in- 
dependent observational probes. This implies that there 
is no strong evidence rejecting the presence of the > 30- 
Mpc ionized bubble at z = 7. 

6. CONCLUSIONS 

We have identified 22 z-dropout galaxy candidates in 
the 0.4 deg2 area of SDF and GOODS-N down to y = 26 
with deep ~ 30-hour Subaru y-band images and Subaru 
and HST legacy imaging data. One out of 22 z-dropout 
galaxies in the SDF has a spectroscopic redshift of z = 
6.96 determined from Lya emission. We have derived 
the bright-end UV LF of galaxies at z = 7. Based on 
our bright z-dropout galaxies as well as faint z-dropout 
galaxies obtained by the recent HST/WFC3 studies, we 
have constrained the early stage of galaxy formation and 
photon budget of cosmic reionization at z = 7. We have 
also discussed the distributions of our z-dropout galaxies. 
The major results of our study are summarized below. 

1 . We find that our bright-end UV LF shows a decrease 
from z 6 to 7. This decrease cannot be explained 
by statistical errors or cosmic variance that is less than 
a factor of < 2 in number density for a 0.2 deg^ area. 
The best-fit Schechter parameters of z = 7 galaxies are 
= 0.69^0^5 X lO^^Mpc^^^ M^y = -20.10 ± 0.76 

SFR thus obtained is 19.7 AIq yr^""^ for the Lya emitting 
dropout. 



mag, and a — —1.72 ± 0.65. Our Schechter parameter fit 
results reject no evolution of UV LF from z = 6 to 7 at 
the > 95% confidence level. A more dominant decrease 
of L* than 0* is preferable from z = 5 — 6 to 7 in the 
error contours of Schechter parameters. 

2. The cosmic SFRD drops from z = 2 — 3toz = 7by 
a factor of at least > 6. It is likely that the cosmic SFRD 
decreases roughly by a factor of ^ 10, but not larger than 
~ 100. This suppression of cosmic SFRD implies that we 
would be witnessing the early phase of galaxy formation 
history. 

3. If we assume the properties of low-z star-forming 
galaxies including the spectral shape and escape fraction 
(/osc — 0.05), we find that the ionized photon produc- 
tion rate of z = 7 galaxies falls below the hydrogen IGM 
recombinatio n rate predicted by the analytic models of 
iMadau et al.l ()1999[ ) even in the homogeneous Universe 
(C'hii = 1) at the ~ 95% {2a) confidence level. Although 
it implies that the Universe cannot be totally ionized by 
only galaxies at z = 7, but we think that properties of 
galaxies at z = 7 are just different from those at low red- 
shifts with, e.g., a larger escape fraction (/esc ^ 0.2), a 
lower metallicity, a flatter IMF, and/or less dust extinc- 
tion. In either case, our observational constraints imply 
that the Universe at z ~ 7 is close to being in balance 
between rates of ionizing photon production and IGM 
hydrogen recombination. 

4. z-dropout galaxies may be strongly clustered both 
in SDF and GOODS-N. We find that the distribution of 
z-dropout galaxies in SDF appears to be a filamentary 
shape which extends up to 60 Mpc in projection, and that 
the z = 6.96 dropout galaxy with a Lya line is located 
at the center of the overdense region consisting of the 4 
UV brightest dropout galaxy candidates. This implies 
that there may exist a well-established ionized bubble 
made by the 4 UV brightest dropout galaxies, and that 
the ionized bubble might help to transmit the Lya line 
in IGM at z = 7. 
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TABLE 1 

Summary of Imaging Observations and Data 



Field 


Band 


CCD 




PSF size'' 


Area 




Date of Observations and Note 








(sec) 


(arcsec) 


(arcmin-^) 


(3(T AB mag) 




SDF 


y 


MIT-LL 


50614 








2007 Apr. 17-20, (2003 Jul. 2, 2004 Mar. 19-20)=^ 




y 


Hamamatsu 


43901 








2009 Apr. 24-26 




y 


(total) 


94515 


0.79 (0.99) 


810 


26.4 






z' 


MIT-LL 


104069 


0.91 (0.99) 


810 


27.7 


Data from Morokuma et al. in preparation"^ 


GOODS-N 


y 


MIT-LL 


89308 








2006 Apr. 1, 3, 2007 Apr. 17-18,20 




y 


Hamamatsu 


29101 








2009 Apr. 24-25, 27 




y 


(total) 


118409 


0.87 (0.87) 


758 


26.2 






z' 


MIT-LL 


39150 


0.87 (0.87) 


758 


26.9 


2007 Apr. 19, (2001 Apr. 20, 22, 
















2003 Apr. 6, 2004 Mar. 16, 2006 Fob. 23-24)== 



^ The FWHM of PSFs in the reduced image. The values in parenthesis indicate the FWHM of PSFs that arc matched with broad-band images in 
each field. 



^ The limiting magnitude defined by a Z(J sky noise level in a l" .8-diameter circular aperture. 
^ The observation dates in parenthesis are those of the Subaru archival data taken by the other teams. 
See Jl 
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TABLE 2 
Limiting Magnitudes 



Band GOODS-N SDF 





(1) 


(2) 


If 


26. 9I 






27.5 


28.9 


^435 


28.4 




Y 


26.2"*^ 


28.1 


Veoe 


28.6 




R 


26. 


28.4 


I or i' 


26.5 


28.1 


1775 


27.9 




2850 


27.6 






26.9 


27.7 




26.2 


26.4 




26.1 




"14.5 


25.5 




ms.s 


23.5 




ms.o 


23.4 




m24 


21.2 





Note. — Three sigma limit- 
ing magnitudes in GOODS-N (1) 
and SDF (2). The magnitudes 
are defined with a l".8-diameter 
aperture for the ground- based 
optical images {UBV Rizy), ex- 
cept for the GOODS-N U, V, 
and R. data. Wc apply a 0".2- 
diamctcr aperture for the HST 
images (B435, Veoe, '^775, ^sso), 
and a 3" .0-diamcter for the 
Spitzer 3.6/im, 4.5/^m, 5.S/J.m, 
8.0/J-m, and 24/j.m images (ms.e, 
^4.5, "i5.8, ms.o, "124)- The 
hmiting magnitudes of Spitzer 
data include the offsets of aper- 
ture corrections (see text). 
^ The detection limits of 
ground-based GOODS-N U, V, 
and R, images are defined with 
a 3" .O-diamcicr aperture due to 
the relatively poor seeing sizes 
{FWHM ^ 1".5; sec the text). 

The 4ct detection limits of y 
magnitudes are 25.9 (GOODS- 
N) and 26.1 (SDF). The 2a up- 
per limits of z magnitudes are 
27.4 and 28.2 in GOODS-N and 
SDF, respectively. 
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TABLE 3 
2 = 7 Galaxy Candidates 



ID U B V B. I/i z y ?^i3.6 ^4.5 ^ ^ y z^poc Note 



GOODS-N 



GOODSN152505° 


> 


27.3 


> 28.0 
(> 28.4) 


> 26.7 
(> 28.6) 


> 


27.0 


> 27.0 
(> 27.9) 


> 27.4 
(> 27.6) 


25.50 


> 23. 4*= 


> 24.2*' 


> 1.90 


Source conf. in Spitzcr band 


GOODSN206815 


> 


27.3 


> 28.0 


> 26.7 


> 


27.0 


> 27.0 


> 27.4 


25.52 






> 1.88 




GOODSN201340 


> 


27.3 


> 28.0 


> 26.7 


> 


27.0 


> 26. 9*= 


> 27.4 


25.69 






> 1.71 


Source confusion 


GOODSN104059 


> 


27.3 


> 28.0 


> 26.7 


> 


27.0 


> 27.0 


> 27.4 


25.79 






> 1.61 




GOODSN108036° 


> 


27.3 


> 28.0 
(> 28.4) 


> 26.7 
(> 28.6) 


> 


27.0 


> 27.0 
(> 27.9) 


> 27.4 
(> 27.6) 


25.82 


> 24.8'' 


> 24.4'' 


> 1.58 


Source conf. in Spitzcr band 


GOODSN203550 


> 


27.3 


> 28.0 


> 26.7 


> 


27.0 


> 27.0 


> 27.4 


25.83 






> 1.57 


Spurious source? 


GOODSN134896 


> 


27.3 


> 28.0 


> 26.7 


> 


27.0 


> 27.0 


> 27.4 


25.86 






> 1.54 




SDF 


SDF63544 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.02 


25.42 






1.60 6.96 


Spectrum in live et al.l 120061 


SDF83878 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.04 


25.49 






1.55 




SDF46975 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.48 


25.51 






1.97 




SDF76507 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.11 


25.59 






1.52 




SDF123919 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.51 


25.79 






1.72 




SDF77202 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.50 


25.86 






1.64 




SDF75298 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.50 


25.87 






1.63 




SDF20911 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.41 


25.89 






1.52 




SDF121488 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.49 


25.97 






1.52 




SDF84539 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


28.16 


25.98 






2.18 




SDF16416 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


> 28.2 


25.99 






> 2.21 




SDF64206 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


> 28.2 


26.07 






> 2.13 




SDF107344 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


> 28.2 


26.08 






> 2.12 




SDF136726 






> 29.4 


> 28.6 


> 


28.8 


> 28.6 


27.59 


26.08 






1.51 




SDF41484 






> 29.4 


> 28.6 


> 


28.8 


> 28. 2*= 


28.05 


26.09 






1.96 


Source confusion 



Note. — The upper limits in the ground-based images, UBV Rlizy, are defined by the 2(T level. The numbers in parenthesis are 3cr upper limits of 
HST/ACS -B435, Vgoq, 2775, and 2350 bands. 

^ Our candidates falling in the central ~ 160 arcmin^ field of GOODS-N with HST and Spitzer images. There are no obvious counterparts in IRAC 5.8/im, 
S.O^im, and MIPS 24/im bands (see Figure[^. The 3(7 upper limits in these bands are (nis.g , 'Tig.o , ^24) — (23.5,23.4,21.2). 
^ Fluxes of these objects arc contaminated by close bright objects on the sky. 
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TABLE 4 
UV Luminosity Function at z = 7 



0* 

(10-3Mpc-3) 
(1) 


(mag) 
(2) 


a 
(3) 


Mag. Range 
(mag) 
(4) 


„obs 

(10-»Mpc-3) 
(5) 


Puv 

(lO^Berg s-i Hz-iMpc-3) 
(6) 


^uppcr 

(lO^Berg s-^' Hz-iMpc-3) 
(7) 




-20.10 ±0.76 


-1.72 ±0.65 


-22.3 < M < -18.0 


cy (j + 2.3 

^■^-1.3 







Note. — (l)-(3): Best-fit Schcchtcr parameters. The values of 0* and are given in units of 10 '^Mpc ^ and AB magnitude, respeetively. The 

reduced o^' the fitting is 0.13. (4): Magnitude range of UV LPs that are used for the fitting. (5)-(6): Number density (in lO^'^Mpc""^) and UV luminosity 
density (in 10^'^erg Hz~^ Mpc"""^) calculated with the best-fit Schechter parameters down to the limit of UV magnitude, Muv ^ defined by the 

HST/WFC3 studies. (7): The upper limit UV luminosity density which is the integral of the best-fit Schechter function down to Muv = oo. 



